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We show that a consistent hydrodynamic description of soft-hadronic 
one- and two-particle observables (the HBT radii) studied in the relativistic 
heavy-ion collisions at RHIC may be obtained if one uses the Gaussian 
energy density profile as the initial condition. The transverse-momentum 
spectra, the elliptic flow coefficient i>2, and the pionic azimuthally sensitive 
HBT radii are successfully reproduced, which hints that the long standing 
HBT puzzle has been solved. 

PACS numbers: 25.75.-q, 25.75.Dw, 25.75.Ld 

The intermediate stages of relativistic heavy-ion collisions are commonly 
described in the framework of the hydrodynamic models [H [H El IH El 
|7j- Such models turned out to be very successful in reproducing the one- 
particle measurements such as the transverse-momentum spectra and the 
elliptic flow coefficient V2- In this context, the early starting time of the 
hydrodynamic evolution that was needed to describe the data as well as the 
explanation of the quite large experimental values of V2 were interpreted 
as the evidence that the quark-gluon plasma created at RHIC is a strongly 
interacting system [SI. On the other hand, the hydrodynamic models failed 
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Fig. 1. The transverse-momentum spectra of pions, kaons and protons for the cen- 
trality bin c = 0-5% (upper panel), c — 20-30% (middle panel), and the elliptic flow 
coefficient for c = 20-40% (lower panel), plotted as functions of the transverse 
momentum and compared to the RHIC Au+ Au data [H [10] . 



in reproducing the experimental two-particle measurements, such as the 
pion correlation functions. The most distinct example is the ratio of the so 
called HBT radii R Q ut an d -Rsidc ~~ in the typical hydrodynamic calculations 
it comes out too large, exceeding the experimentally measured values by 
about 20-50%. 

In our recent work [TT], we have found that the consistent description 
of one- and two-particle observables in the soft region is achieved within 
the hydrodynamic model if one modifies the initial conditions - the initial 
energy profile obtained in most cases from the optical Glauber model should 
be replaced by the Gaussian profile. Other ingredients of our model, such as 
the use of the semi-hard equation of state, are also important for the overall 
good agreement, however, taken into account alone they are not sufficient to 
get the satisfactory results. In this aspect, our findings are complementary 
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to the observations made recently by Pratt (Ref. [12] and these proceedings 

Our approach is based on the 2+1 boost-invariant inviscid hydrody- 
namics |14l [T5l [16] followed by the statistical-hadronization model THER- 
MINATOR [T7]. The initial energy-density profiles are obtained from the 
Monte-Carlo Glauber model GLISSANDO [18], which includes the eccen- 
tricity fluctuations [El [201 ED E21 [23]. The results of those calculations 
are approximated by the Gaussians. The thermal event generator THER- 
MINATOR simulates hadron emission from the freeze-out hypersurface de- 
termined by the hydrodynamic calculation. The single freeze-out scenario 
[241 [25] is assumed with the universal final temperature Tf =145 MeV0- 
Besides the final temperature our model has essentially two additional pa- 
rameters: the initial temperature Tj, fixing the absolute normalization, and 
the initial time for the start of hydrodynamics, tq = 0.25 fm. For each 
centrality class, the widths of the two-dimensional initial Gaussian energy 
distribution are obtained directly from the GLISSANDO simulations. 

Figure Q] shows our results for the transverse-momentum spectra of pions, 
kaons, and protons for the centrality classes c = 0-5% and c = 20-30%. 
Fig. [1] presents also our results for V2 for the centrality c = 20-40% , plotted 
as functions of the transverse momentum. The spectra and the elliptic 
flow are compared to the RHIC data [9l [10]. We observe a very good 
agreement between the model predictions and the data. The small excess 
of the theoretical proton vi above the data may be attributed to the lack of 
rescattering in the final state. 

Figure [2] presents our results for the pion HBT radii i? s ide , -Rout, -Rkmg, 
and the ratio Rout/Rside for central collisions, compared to the RHIC data 
[27] . The left panel shows our best results obtained with the traditional 
Glauber initial condition [28J, while the right panel shows the results ob- 
tained with the Gaussian initial condition [11]. One observes that a very 
good agreement between the data and the theoretical model predictions is 
achieved in the case where the Gaussian initial condition is used. We note 
that the calculation of the radii does not introduce any extra parameters. 
All the characteristics of the emitting source were already fixed by the fits to 
the spectra and V2- Finally, in Fig. [3]the results describing the azimuthal de- 
pendence of the HBT radii are plotted [29]. Here R 2 ^) = R^ + 2Rj cos (2^). 
Again, we observe a very good agreement between the data and our model 
for different centralities and different average momenta of the pion pairs kx- 

Our results indicate that it is possible to achieve a consistent hydrody- 
namic description of one- and two-particle soft hadronic data collected in the 

1 The relatively high decoupling temperature may be connected with the presence of 
the instabilities at the phase transition connected with the diverging bulk viscosity 
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Fig. 2. The pion HBT radii i? S idc , -Rout , -Riong, and the ratio i? ou t/-Rsidc for 
central collisions, shown as the functions of the average momentum of the pair and 
compared to the RHIC Au+Au data [27]. Left: our best results obtained with the 
Glauber initial condition |28j . Right: our results obtained with the Gaussian initial 
conditions (i.e., with the Gaussian approximation to the the Glauber simulations) 

na. 



RHIC heavy-ion experiments at the highest beam energies. In particular, 
it is possible to describe simultaneously the transverse-momentum spectra, 
the elliptic flow coefficient 1)2, and the HBT correlation radii (the prelim- 
inary results show that the correlations of the non-identical particles are 
also well reproduced in our model). Our finding suggests that there exists a 
solution to the long standing RHIC HBT puzzle understood as the failure of 
the relativistic hydrodynamics in reproducing consistently the spectra and 
the HBT radii. 

The Gaussian initial profile leads to a faster development of the initial 
transverse flow which makes the system evolution shorter. At the same time 
the transverse size of the system at freeze-out is slightly larger (as compared 
to the standard Glauber scenario). These two effects put together lead to 
the desired reduction of the ratio i? out /-R s idc • In the approach discussed in 
Refs. |12tll3j. the crucial role of the faster initial acceleration |144I30| is also 
emphasized. However, in [13] this effect comes from the very early start of 
the hydrodynamics, tq = 0.1, fm and from the inclusion of viscous terms 
which increase the transverse acceleration with respect to the longitudinal 
one. We stress, however, that our approach remains the only one which 
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Fig. 3. Results for the RHIC HBT radii and their azimuthal oscillations. For each 
value of -ZV part on the horizontal axis the experimental points (filled symbols) and the 
model results (empty symbols) are plotted. The points from top to bottom at each 
plot correspond to hp contained in the bins 0.15-0.25 GeV (circles), 0.25-0.35 GeV 
(squares), and 0.35-0.6 GeV(triangles). The top panels show R% ut , Rg idcQ , and 
i?f ong Q, the bottom panels show the magnitude of the allowed oscillations divided 
conventionally by R^ idc • 



describes correctly the HBT radii and the elliptic flow, as the studies of 
Refs. \12\ [T3] are restricted to the cylindrically symmetric case. 

Clearly, the issue of the relative importance of the initial energy profile, 
the initial pre- equilibrium flow, and the viscous effects should be addressed 
in further investigations. The recent calculation [31] shows that the pre- 
equilibrium flow is inevitable as it comes from the energy-momentum conser- 
vation laws. On the other hand, the specific amount of the pre-equilibrium 
flow required to correctly describe the data depends on the initial energy- 
density profile. In [32] we have assumed that the start of the hydrodynamic 
evolution may be delayed to r = 1 fm. The earlier evolution consists of the 
free-streaming (in the proper time interval: 0.25 fm < r < 1 fm) followed 
by the sudden equilibration and transition to the hydrodynamic regime (at 
r = 1 fm). We have used again the Gaussian initial conditions and shown 
that they lead the sufficient amount of the pre-equilibrium transverse flow 
at t = 1 fm, when the hydrodynamic evolution starts, i.e., the FS+SE ap- 
proach of Ref. [32] describes the data equally well as our standard approach. 
This behavior suggest some interesting universality of the Gaussian initial 
conditions. Of course, the open question remains to find the microscopic 
mechanism leading to such a form of the initial conditions. 
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